class: center, middle, inverse, title-slide # Übung 2: Komplexe Vektordaten ## CAS FAB: Räumliche Daten in R ### Nils Ratnaweera ### Forschungsgruppe Geoinformatik ### 2021-11-30 --- ## Vorbereitung <!--1--> Erstelle ein neues R Script mit dem Namen `Uebung_2.R` und lade darin die libraries `sf` und `tmap` ```r library(sf) library(tmap) ``` --- ## Übung 2.1 <!--10--> - Suche dir die Gemeindegrenzen der Schweiz. - Drei nützliche Adressen hierfür sind: - [opendata.swiss](https://opendata.swiss/) - [map.geo.admin.ch](https://map.geo.admin.ch) - [swisstopo.admin.ch](https://www.swisstopo.admin.ch) - Wenn du die Wahl hast, versuche das File als Geopackage herunterzuladen. Ansonsten als Shapefile oder als File Geodatabase - Entzippe das File (sofern nötig) und schau dir den Inhalt an -- ## Lösung https://data.geo.admin.ch/ch.swisstopo.swissboundaries3d-gemeinde-flaeche.fill/shp/2056/ch.swisstopo.swissboundaries3d-gemeinde-flaeche.fill.zip Shortlink (für diesen Kurs): https://bit.ly/3CaAj5W --- ## Übung 2.2 Importiere das den Datensatz mit `read_sf()` und speichere den output in der Variabel `gemeindegrenzen` -- ## Lösung ```r gemeindegrenzen <- read_sf("_data/original/ch.swisstopo.swissboundaries3d-gemeinde-flaeche.fill/swissBOUNDARIES3D_1_3_TLM_BEZIRKSGEBIET.shp") ``` --- ## Übung 2.3 <!--5--> Betrachte den importierten Datensatz in der Konsole. Was für Informationen kannst du entnehmen? -- ## Lösung - Die Koordinaten sind im Bezugssystem `CH1903+ LV95` abgespeichert - Es handelt sich um 159 Gemeinden und 21 Attribute, die meisten unbrauchbar --- ## Übung 2.4 <!--3--> Entferne alle Spalten bis auf `NAME`, `EINWOHNERZ` und `geometry`. -- ## Lösung ```r gemeindegrenzen <- gemeindegrenzen[, c("NAME","EINWOHNERZ","geometry")] ``` --- ## Übung 2.5 <!--4--> Visualisiere die Gemeindegrenzen mit `plot()` und `tmap`. -- ## Lösung ```r tm_shape(gemeindegrenzen) + tm_polygons() ```